Theory of correlated electron transport and inelastic tunneling spectroscopy 
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For a non-superconducting system, the electronic tunneling current through an insulating barrier 
is calculated, including interaction effects. The exact Hamiltonian of the full system is projected 
onto the subspaces of the "left" and "right" leads. In the weak tunneling limit the well-known 
tunneling Hamiltonian is recovered, along with an additional term. This additional term originates 
from the projection of the electron-electron interaction onto each subsystem and corresponds to 
correlated tunneling. It is shown that the tunneling current is determined by — in addition to the 
single-particle density of states — the spin-spin and density-density susceptibilities. The signatures 
of which have recently been observed in several experiments. 

PACS numbers: 73.40.Gk,72.10.-d 



I. INTRODUCTION 

Tunneling experiments were one of the first major con- 
firmations of the Bardccn-Coopcr-SchriefFcr (BCS) the- 
ory of superconductivity.^ These experiments showed a 
dramatic suppression in the conductance near the Fermi 
energy. The theoretical explanation of this was the for- 
mation of a gap in the quasi-particle density of states of 
the superconductor.-^ This essentially began the relation- 
ship linking a tunneling current in a many-body system 
to the single-particle density of states. 

To describe tunneling at a general many-body level, 
Cohen et al? later introduced the tunneling or transfer 
Hamiltonian, consisting of two disjoint, fully interacting 
Hamiltonians connected by a tunneling term. This tun- 
neling term removes electrons from one system and cre- 
ates them in the other. Although phenomenological, the 
tunneling Hamiltonian approach has become and remains 
the de facto standard for theoretically interpreting tun- 
neling experiments in condensed matter, across a wide 
range of systems. Several attempts have been made to 
put many-body tunneling, with or without the tunneling 
Hamiltonian, on more rigorous theoretical ground,"' but 
owing to the inherent difficulties of this problem, it still 
remains open. 

Nevertheless, the standard tunneling formalism has 
been used extensively to account for an extraordinarily 
wide variety of experiments. Although, with the ever 
increasing refinement of experimental techniques, exper- 
iments are starting to probe physics not captured by this 
formalism. One area where modifications are needed is 
in the regime of inelastic tunneling, where the tunnel- 
ing electron loses energy through interactions with the 
barrier or other additional degrees of freedom of the sys- 
tem being probed. These effects are typically accounted 
for by simply adding additional transfer terms. '^ Inelas- 
tic electron tunneling spectroscopy (lETS), most notably 
scanning tunneling microscopy lETS (STM-IETS), has 
be used to probe the internal vibrational modes of molec- 
ular adsorbates on surfaces. More recent experiments®"^ 
on single magnetic atoms, finite spin chains, and mag- 
netic substrates have shown evidence of spin related in- 



elastic tunneling. Here, especially for the adatom chains, 
the current-voltage (TV) curves show features where one 
would expect spin excitations of the system, such as the 
singlet-triplet excitation of neighboring pairs. Explaining 
such characteristics is outside the scope of the standard 
tunneling approach, which relates an I-V curve, or con- 
ductance, to single-particle properties, namely the local 
single-particle density of states. Spin excitations, such as 
the singlet-triplet transition, arc manifestly two-particle 
properties, related to the spin (density) susceptibility or 
spin-spin correlation function. Several authors have been 
able to account for the experimental results of these ex- 
periments remarkably well,^'^^ by introducing an inter- 
action effect between the STM tip and adatoms, leading 
to so-called spin mediated tunneling. The microscopic 
origin of such a term and how it should be generalized, 
to say other spin excitations such as magons,^ isn't im- 
mediately clear though. In principle these and other in- 
teraction effects should occur in virtually all tunneling 
experiments.^^ One would like a theoretical foundation 
that is first principles driven such that it could be applied 
to a wide variety of systems and also be generalizable to 
others in a systematic manner. The main goal of this ar- 
ticle is to give such a derivation of the tunneling current, 
in the spirit of the transfer Hamiltonian, that includes 
these interaction or two-body effects. 

In the next section we start with a general interact- 
ing Hamiltonian for a system with a tunneling barrier. 
The original Hamiltonian is then projected onto the low- 
energy (below the height of the barrier) "left" and "right" 
subspaces. This projection is conjectured to capture the 
most salient tunneling processes. From which, the stan- 
dard tunneling Hamiltonian is recovered along with an 
additional term, which describes correlated or interac- 
tion meditated tunneling, stemming from the original 
electron-electron interaction. In Sec. Ill we calculate the 
steady-state tunneling current through the barrier with 
respect to this effective tunneling Hamiltonian. The re- 
sult gives an additional contribution to the current from 
two-particle susceptibilities, as well as the well-known 
contribution determined by the single-particle density of 
states. 



II. DERIVATION OF EFFECTIVE TUNNELING 
HAMILTONIAN 

Here we consider a one-dimensional system for nota- 
tional convenience; the generalization to higher dimen- 
sions is straightforward. For an interacting system in 
the presence of a tunneling barrier, the grand-canonical 
Hamiltonian of the entire system is given in second quan- 
tization by (setting fi, = 1) 



H = Y, dx^ii 



2m 



fi + V{x) ^^{x) 



a, a' 

(1) 

where U{x,x') is the electron-electron interaction (as- 
sumed to be symmetric), /x is the chemical potential, and 
V{x) is taken to be the tunneling barrier, but could also 
contain disorder, lattice fields or other single-particle po- 
tentials, sec Fig. 1. The field operators ^^(a;) and '^^{x) 
obey the standard fermionic anticommuation relations 
{^,{x), ^^.{x')} = 0, {^,ix),^l,ix')} = 5,,„.5{x ~ x'), 
and can be expressed in terms of mode creation and an- 
nihilation operators by 



k 



fccr' 



(2) 
(3) 



where {ak^,a^,^,} = 0, {afc<,,4'^'} = S'y,<T'Sk,k' , and 
ip^^{x) are the exact single-particle eigenstates of (1). 
Next, we consider two related Hamiltonians 



Hl = Y. dx^ljx) 



2m 



^i + Vl^ix) ^i^Jx) 



a, a' 

(4) 
and 



V^ 



Hr = J2 dx ^l^x) - ^ - A^ + V^{x) ^^^„{x) 



^J2l'i^dx'^ljx)^l^,{x')U{x,x')^^,,,{x')^^Jx 



where 



Vdx) 



V{x), for X < a 
V{a), for X > a , 



(5) 



(6) 



V(x) 

n 



t 




FIG. 1. Schematic of a tunneling barrier for the full Hamilto- 
nian, Eq. (1), having a characteristic width d and height Vo- 



V,(x) 




and 



FIG. 2. Potential barrier for Hi,, defined in (6). 



^''^''' ~ \V{a), ioTX<a. 



(7) 



The parameter a is an arbitrary point within the barrier, 
such that V{x) = VL{x)Q{a — x) + Vb,{x)Q{x — a), where 
Q{x) is the Heavyside step function; see Figs. 2 and 3. 
The field operators in (4) and (5) are given by 



k 



R, fccr J 



(8) 
(9) 



where tph ka ^'^'^ V'r ka ^^^ ^^^ exact and complete single- 
particle eigenstates of Hi, and Hr respectively. 

The Hilbert space of H , % can be decomposed into the 
Hilbcrt spaces of Hi, and i^R, 



n = HL + nn + sn, 



(10) 



where SH allows for states outside of either Hi, or Hr, 
e.g. resonant states within the barrier itself. Such possi- 
bilities will not be considered here. One can then define 







FIG. 3. Potential barrier for J/r, defined in (7). 
projection operators 

fl. = X! IV'L.q)(V'L,q 
a 

Pr = X! IV'R,q)(V'R,q|, 



(11) 



which project the original Hilbert space on to the left and 
right subspaces respectively. We assume that the sum of 
projection operators, to a good approximation, forms a 
resolution of the identity in "H, i.e. 

^L + ^R^a- (12) 

The electron field operator of (1) can therefore be written 
as 

^^i^{x) + ^r{x)+5^{x), (13) 

where S'i'{x) is of order of the tunneling amplitude, i.e., 
for states below the barrier height 6'i>{x) ~ 0{Ta.a' ~ 
{^L,a\V\ipR^a'))- (5^(.t) will be neglected from here on, 
as it will ultimately correspond to higher order terms, in 
the tunneling amplitude, in a calculation of the current. 
Projecting the field operator of the original Hamiltonian 
onto the left and right subspaces 

k 

*[,.(^) = pI^U^)Pl = E ^fe.(^)«Lk.: (14) 



and 



k 

^lA^) - pI'^U^^Pr = E ^^-(^)<k.. (15) 



where 



a 



L,k(T 



Idx'lp^^kai^)'^^^)^ 4,k<T 



'^R ko- ~ /'^^V'r kiji^)'^a{^)- From which it readily fol- 
lows that {aL^ka'«U'<T'} = {«R.k^'"R,k'<T'} ^ '5^,<T'4,fc', 

l'^L,k<T''*L,k'o-'l — l'^R,kCT''^R,k'cr'l ~ I'^L^ko-' ^R^kV I ~ 

0, and {aL.k^'aR,k'<T'} = ('/'L.fc^l Vfe'^^')- 1^ S'^^- 
eral the left and right states are not orthogonal; 
(V'l.qIV'r Q') 7^ 0. For states with energy below the 
barrier height (V'L^alV'R.a') ~ 0{Ta,a')- 

As the exact single-particle eigenstates of (1), (4), and 
(5) are assumed to form a complete basics, we can expand 
the original eigenfunctions in terms of the left and right 
basis 



-ipkaix) = E'^fc,fc''^L,fcv(a;) 
fe' 

tpka{x) =^cfk'i'R,k'a{x), 



(16) 
(17) 



Idxipl^k„{x)'iUx), Oj^^k^ = Jdx'tlj^,^^{x)^^{x), and 



where c]:,^, = {iphM^-'k') and cff., = (V'R,fc|'0fe')- The 
expansion coefficients can be obtained from assuming 
an explicit form of the tunneling barrier** or within 
WKB (Wentzel-Kramers-Brillouin) theory for an arbi- 
trary potential. ^■^ For a wide tall barrier and for states 
with energy much less than the barrier height c^ ^, oc 

4,fc' + 0(Tfe,feO and cf-^, oc SkM' + 0{Tk.k')- Thus, apart 
from an overall proportionality constant, to leading order 
in the tunneling amplitude 

^L,a(x)^*L,a(a:;) + O(Tfc^fc0 (18) 

^kA^) - -^rA^) + 0{Tk,k')- (19) 

Therefore, the original fields operators, to leading order 
in the tunneling, can be decomposed into the left and 
right subsystems by 

^„{x) ~ ^LA^) + ^RA^) + o{Tk^k') 

^Ux) ^ *L(2^) + *L(^) + o(rfc,fcO. (20) 

It should be noted that, at this point, because we have 
restricted ourselves to states below the barrier, the field 
operators "ij^ix) and ^j^(x) in Eq. (20), no longer create 
(annihilate) electrons in localized delta-function states, 
but instead in some broaden approximation. We will not 
deal with this technicality and assume that for a high 
barrier, compared to the Fermi energy, to a good approx- 
imation the states below the barrier form a sufficiently 
complete set. 

Expressing the exact Hamiltonian (1) in terms of the 
approximate projected operators (20), and for now ne- 
glecting interactions, gives 

^ ~ J2^^ka - M)4,k<T«L.ka + E(^fe'^ ^ M)aR,ka«R,kcT 

k,cr k.a 

+ E [^fc.fc'^lkaaR.kv + H-C-] 

+ E [^fe,fe','^('=fcV"M)aLkaaR,kV + H-C-]' (21) 

k.k' ,(T 



where terms such as / dx '^l^{x)V {x)'i' j^{x) that only lead 
to a change in the density and energies near the bar- 
rier have been neglected, and ^k.k' — (V'l fclV'R fc')- ^^ 
this level, the fact that the eigenstates of the left and 
right sides are not orthogonal simply leads to a renor- 
malization of the tunneling matrix elements, Tk^k' ^ 
Tk,k> + £,k,k' (cfe' — m)- Even in the approximation S,k.k' ^• 
0, the important and dominate physics is stilled cap- 
tured. Thus, neglecting these terms or equivalently set- 
ting ^fc.fc' = 0, one can still expect to obtain quantita- 
tively correct results, of course this approximation can be 
relaxed. ^'^ Also, this simplification exactly reduces (21) 
to the standard tunneling Hamiltonian. In addition, this 
greatly simplifies calculations of expectation values, i.e 
Green's functions, as now {"^^ ^(a;), ^r ^(a;)} = 0, which 
we will assume from here on. With such an approxi- 
mation we have essentially arrived back at the starting 
point for the standard tunneling Hamiltonian. One may 
wonder why all of the previous formalities were neces- 
sary. For one, it demonstrates what approximations and 



assumptions are made when one is using the standard 
tunneling Hamiltonian formalism, and secondly it will 
allow us to treat the effects of interactions on tunneling 
quantitatively, which we now turn to. 

Normally, interactions are introduced simply by adding 
the appropriate interactions terms to the left and right 
subsystems. This neglects the cross terms generated by 
expressing the original operators in terms of the approxi- 
mate left and right states. As we will see, it is these cross 
terms that are responsible for the electronic inelastic tun- 
neling properties. Expressing the interaction term of (1) 
in terms of the projected operators, Eqs. (20), leads to 
(see the Appendix A for details) 



H « _ff L + Hj{ + Hi 



(22) 



where Hi^, Hji are the fully interacting Hamiltonians of 
the left and right subsystems 



H^^Y. dx^lj 



2m 



^lJ^) + o E / dxdx'M'ljx)¥^,^,{x')U{x,x')^'^^,,{x')^'^^,{x), 



Hn = Y. dx^l^ix) 



2m 



*R,.(-^) + o E / dxdx'^ljx)^l,,ix')U{x,x')^^^„,{x')^^,Jx), (23) 



and the tunneling or transfer part is given by 



i/T = E / '^^ [T{x)'il,ix)^n,A^) + H.C.] +Y,J dxdx' U{x, x')[^l^{x)^^^,{x) [n^{x') + n^{x')\ + H.c}, 



(24) 



where n{x) = I]^. *J.(a:)*a(a;). The first term of (24) 
is the standard single-particle tunneling, which comes 
from the kinetic energy term of the original Hamiltonian. 
The last term is from the interaction of the original sys- 
tem and is normally not considered in most treatments 
of tunneling. This term involves a interaction mediated 
transfer of electrons, or correlated hopping in a lattice 
model. ^"^ Although, this term might be small, as it in- 
volves the correlated transfer electrons, it is ultimately 
responsible for the inelastic tunneling properties. 



III. EVALUATION OF THE TUNNELING 
CURRENT 

In the presence of an applied electric field a net current 
will fiow through the barrier. Here we will calculate the 



I 

steady-state current in the limit of weak tunneling be- 
tween the left and right sub-systems. As usual, ^^ we will 
assume that each subsystem is separately in thermody- 
namic equilibrium, where the chemical potential of each 
side differs only by the applied voltage eV . The drop in 
the chemical potential is assumed to occur entirely within 
the barrier. Because the total particle number commutes 
with the effective Hamiltonian, the current is defined to 
be proportional to the expectation value of the time-rate- 
change of the number of electrons in the left (or right) 
sub-system. By Heisenberg's equation of motion (e > 0) 
the current operator is given by 

/ = -ed^Ni^ = -ie [iVL, H] = -ie [TVl, -ffx] , (25) 

where TVl = ^a I dx'^\^^{x)'^j^^{x). Evaluating the 
commutator gives 



/ = — ze 



; J dx [T{x)^l^{x)^^^A^) - H.C.] - *e ^ / dxdx' U{x, x'){^ljx)^^^A^) [Mx') + iMx')] - H.c}. 



The nonequilibrium expectation value of (26) gives the 
measured tunneling current. Here we obtain this expec- 
tation within linear response theory, treating the tunnel- 
ing term, Hj^, as the perturbation. To leading order in 
the tunneling, the current is given by the thermal expec- 
tation of 



(27) 



/-* / dt'{[HT{t'),m])^^^, 



where 0(t) = e'"°*Oe~'"'>\ Hq = Hi^ + H^, and we 
have assumed that as i — > — oo the two sub-systems are 
completely decoupled, i.e. the height of the barrier is in- 
finitely large. 

Evaluating (27) for a non-superconducting system and 
assuming local tunneling, i.e., the wave functions of one 
side are spatially localized about a point x, as is the case 
for an STM, the total current can be written as (See 

I 



(26) 
I 

Appendix B for a complete derivation.) 

I^h + h- (28) 

The first term gives the standard expression for the tun- 
neling current. 

/i = 27re|rp^ / duj pi^,^{x,uj + eV)pj:i„{x,uj) 

X [nF(w)-nF(w + el/)], (29) 

which relates the single-particle local density of states 



p^(a;,cj) = IniG;'=*(x,a;), 



(30) 



to the total current, where G"'°* is the retarded single- 
particle Green's function and n-piyi) — (exp(/3aj) 4- 1)~^. 
The second term is related to the electronic inelastic tun- 
neling and is in general given by 



/2 = 

OO 

dJ^Y. J du;du;'p^Jx,iv + el/)pR_,,(x, w'){[Xf' ,?,?,?' (a;, w' - w - eF) + x^-,^,^,y (x, w' - u;)]np{u;')[l - ripiu; + eV)] 
-heTT^/^^ / dujdJp^^{x,uj + ey)p^,^(x,J)p^l,^^(x,uj + J^eV) 

-CTTt/^^ / dujdJ p^-^{x,uj^ey)p^^^(x,J)p\„;s{x,uj+J) 

X In^iij: \- eV)n^(ui')\Y - n^{uj ^ J)\ - \Y - n^{uj + eV)\ \Y - n^{J)\n^{uj + J)\, (31) 

where X^,,o^,o,,aA^^'^) = {^\,aS^^^)'^ x,o^^^^^)^\,a^^^M^ x,oA^^^)) h^ and p%[x,^) = -ilmG;,i^(x,w) is the 
two-particle density of states, with 

G^l^{x,t) = -i0(t)({vl/,(x,t)vI/-(x,t),vl/;(O)vl/t(O)})^, (32) 

and a is the opposite of g. Note that x is not a time-ordered expectation, but is instead, apart from thermal 
factors, related to a two-particle spectral function. For a system with weak or no correlations on one side and a 
spin-independent featureless density of states near the Fermi energy ep and neglecting the contributions from the 
two-particle density of states. 



/2«e[/VL(eF)pR(eF)5I \ du:du:'^,^„^,^„,(x,uj-J)[n^(Lo){\~n^{J ^ey)\ - n^{LO + ey){\ - n^{J)\^. (33) 



Equation (33) can be rewritten using 

Y,Xa',a..,a'{^.t) = 25r,(x, t) + inR(x,t) 



(T,(T' 



2(sr(2;, t) ■ Sr (x, 0)) + -{n^ix, t)n^{x, 0)), 



(34) 



where {s^{x,t) ■ Sj^(a;, 0)) and {h^{x,t)n^{x^O)) are the spin-spin (density) and density-density correlation functions 
respectively, giving 



I2 = 2eU Pi^i^p) PYii^p) / dujduj' Syi{x,uj — uj')<np{uj)[l — np{uj' + eV)] — np{uj + eV)[l — np{uj')\ > 



eU^p^{eF)p^ieF) 



dwduj' nj^{x, uj - J)\np{bj) [l - np{J + eV)\ - np{uj + eV) [l - np{J)\ \. (35) 



For spin systems, such as those of Refs. [6-8], the den- 
sity fluctuation term can be safely neglected. This re- 
duces Eq. (35) to the form of those proprosed in Refs. 9 
and 10, which have accurately described the experiments 
of Refs. 6 and 7. But in general, Eq. (31) gives the total 
correlated transfer contribution to the tunneling current 
and can be used in a broader range of experiments, which 
will be the subject of a future publications. 

Within the approximations leading to Eq. (35) and in 
the low temperature limit, fcgT <C ep, the differential 
conductance G = 9y/ is 

G{ey) = 47re2|T|VL(eF)pR(eF) + le^V" p^{eF)p^{eF) 

X < doj Sji{x, oj) [np{uj + eV) + np{uj — eV)\ 



+ - / dujIl^{x,Lu)[np{LO + eV) + np(a; - eV)] >. 

(36) 
The second derivative of the current gives 

dl.1 = 2e^U^p^{eF)p^{eF)[s^ix, eV) - S^{x, -eV) 

+ \[n^{x,eV)-U^{x,-eV)\]. (37) 



IV. 



OUTLOOK 



As mentioned in the introduction and demonstrated 
in the body of the article, in principle all tunneling ex- 
periments, at some level, probe two-particle properties 
induced by interactions effects. Although, by their very 
nature these effects can be small, when compared to 
the single-particle terms. It takes a well engineered ex- 
periment, such as those of Refs. [6-8], to separate the 
contributions. Another system where such effects have 
probably been observed, is in transport through a quan- 
tum point contact (QPC). It is believed that interactions 



r 



arc responsible for the "0.7 conductance anomaly" of 
QPCs.^^ A lot of theoretical work has gone into explain- 
ing this effect. For instance, in Ref. [17] a generalized An- 
derson model, which effectively contains correlated trans- 
fer of electrons, has been applied with some success to 
this problem. In the previous sections we have calculated 
the contribution from interactions to the many-body tun- 
neling current within an "STM-like" geometry. A natural 
extension would be to apply the ideas developed here to a 
quantum-dot geometry or QPC. Such a derivation would 
lead to a generalization of the well-known Meir-Wingreen 
expression.^** 

Furthermore, one could use the results presented here 
to design experiments to probe the two-particle prop- 
erties of other important and interesting systems. For 
example, how these properties change as one approaches 
the superconducting instability. This could give informa- 
tion about the pairing mechanism of high-temperature 
superconductors. 
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Appendix A: Details of the derivation of the 
tunneling Hamiltonian 

Here we give the details of the derivation of Eq. (22). 
The kinetic energy term was done in the body of the pa- 
per, here we will deal with the interaction term. Express- 
ing the interaction term of (1) in terms of the projected 
operators, Eqs. (20) and keeping only terms to leading 



order in the tunneling amplitude gives 



^Y. fdxdx'^l{x)¥^,{x')U{x,x')^^,{x')^^{x) = 
^J2 1 dxdx' ^ljx)^l,,{x')U{x,x')^^^,,{x')^^j3 



^J2j d^dx' ^l^M'^i^a' i^')U{x, XO^R,,,, (X')*R^, 



Left interaction 



Right interaction 



Y, j dxdx' ^lA^)^la' {^')U{x, x')*R,.' {x')^i.A^) + E / ^^^^' *U(^)*R,.' i^')U{x. a;')*L,.' {x')^kA' 



Direct intcracti( 



Exchange interaction 



E I dxdx' t/(x, x'){<,(a:)vI/R^.(a;) [nL(x') + nR(x')] + H.c.} +0(r2,,), 



(Al) 



Tunneling 



where h{x) ~ ^^ '^'^{x)'^ „{x) . The exchange and direct 
terms don't contribute to tunneling current, i.e. they 
both commute with the number operator of each side, 
see Sec. III. and thus will be neglected. The full effective 
Hamiltonian then reduces to that given by Eq. (22). 



Appendix B: Details of the derivation of the 
tunneling current 

Here the details are given for the evaluation of the 
tunneling current, Eq. (27). Although the evaluation is 
quite lengthy and tedious, a lot of details are included 
for completeness. 

Defining A^{x) = T {x)^'l,^{x)^ ^,^{x) and 

B„{x,x') = ;7(x,x')*[,,(a;)*R^,(a;)[nL(x')+riL(x')] for 
notational convenance. Note that 

{A,{x))h,^{B,{x,x'))h,=Q 
(i.(x)i,,(x'))Ho = {Mx,x')B^,{x'\x'"))h^^=Q 



and 



{A„{x)B,.{x',x"))h,=Q. 



Evaluating the commutator in Eq. (27), the total current 
can be written as / = /i + /2 + /a , where 



/i =2e / dt'Y / dxdx Re {AI,{x' ,t')A„{x,t)) 

- {A,,{x' ,t')Al{x,t)) 



(Bl) 



t 
/2 = 2e / dt' y I dxidx2dx3dx4 



X Re^{Bl,{xi,X2,t')B^{x3,X4,t)) 

- {B,, ixi,X2,t')Blix3,X4, t))] , (B2) 



and 



t 
I3 ~ 2e / dt' y / dx'dxidx2 



Re[{Al{x',t')B^{xi,X2,t)} - {A^,ix',t')Bl{xuX2,t)) 

+ {Bl{xi,X2,t')A,{x',t)) - {B,,{xuX2,t')Al{x\t)) 

(B3) 

The first two contributions will be evaluated in the fol- 
lowing sections, while the third term, /a, will be ne- 
glected. As it contains two-particle correlations with 
three equal times. It would be expected that these cor- 
relations are subdominant by this restriction of phase- 
space. 

1. h 

To evaluate the expectation values in Eq. (Bl) it is 
useful to introduce the Keldysh contour-ordered single- 
particle Green's function^^ 

G^{x,tX) = -z(Tc,*^,,(x,i)*t,_,(x,i'))^,^., (B4) 

where X £ {L, R}. The quantities of interest for this 
work are the "lesser" (<) and "greater" (>) functions 



given by 

(B5) 
and 

= -z fdoj e-'-(*-*') [1 - np{Lo)] p^Jx, cj), (B6) 



where rip (a;) = (e'^" + 1)~^ is the Fermi distribution with 
inverse temperature /3 — {kBT)~^ and p]^{x,uj) is the 
single-particle local density of states, obtained from the 
imaginary part of the retarded Green's function, 



Pai^,^) 



-lmGl'''ix,Lu). 



(B7) 



In the limit were the tunneling is spatially localized 
about a single point, as it for an STM, where the wave 
functions of the tip are exponentially localized on the 
atomic scale about the location of the tip, the spatial in- 
tegrals in (Bl) can be approximated by the value of the 
functions centered at the tunneling point. This is equiva- 
lent to the, common, assumption that the tunneling ma- 
trix elements Tk,k' = ('/'L,k|2^|V'R,k') have no momentum 
dependance. The needed correlation functions from (Bl) 
are then (assuming spin is conserved during tunneling) 



{Al{x,t')A„ix,t))Ho^ 

(B8) 



ir|2(*, 



and 



{AJx,t')Alix,t))Ho = 

(B9) 



By using Eqs. (B5) and (B6), the well-known result re- 
lating the tunneling current to the single-particle local 
density of states is recovered, 



h = 27re|Tp^ / duj pi^,^{x,uj + eV)pj:i„{x,uj) 
X [n^iuj) - np{uj + eV)] . (BIO) 



h 



The remaining contribution to the tunneling current, 
I2, contains many high order correlations functions, in- 
cluding three particle propagators. The evaluations of 
which becomes quite tedious and lengthy. But as this is 
directly related to this work, we will cover these in some 
detail. 

From Eq. (B2), the needed correlations for I2 are 



{BAxi,X2,t')B^{x3,X4,t))H^ = 

U{XI,X2)U{X3,X4) 



;*L,.(a;i, i')*La(2;3, t))H^ {n^ix2,t')^iAxi,t')^j,Ax3,t)n^ix4, t)) ^^ 
[nLix2,t')^i^Axi,t')^ljx3, t))H^ (*|,,^(;ri, t')*R.^(a;3, t)nAxi. t)) h^ 
;*L..(a;i, i')*La(2^3, t)nA^i, t))H^ (np (X2, i')*L(a^i, i')*R,.(a^3, t)) h^ 



and 

{B^{xi,X2,t')Bl{x3, 2:4, t))u. 



C/(xi, a;2)C/(x3, x^) [(*[^^(xi, t')^L(^2, t')^L(a;4, t)^i^,A^i. t)) h^ (*R,.(a^i, i')*La(a;3, t)) 

+ (*[,.(a;i: i')*L,a(a^3, 0)ffL (*R,.(a;i, i')^R(a^2, i')"R(a;4, t)'^iA^3,t)) h^ 
+ (*[,.(a;i, i')"L(a;2, i')*L..(a;3, i))if^ (*r,,,(.ti, t')^R(a^4, i)*R,.(a;3, t))H^ 
+ (*[..(a;i, i')"L(a;4, i)*L..(2;3, 0)hl (*R,<r(a;i, i')"R(a^2, i')*R.a(a;3, t))H 



(B12) 



Neglecting all two-particle correlations that have three equal times, as these are expected to be sub-dominate from 
the reduction of phase space, leads to 

U{XI,X2)U{X3,X4) \{h^{x2,t')^^^^{xi,t')^ljx3,t)h^{x4, t)) „ (*|^_^(xi , t')*R,^ (a^s, t))H^ 



(B13) 



+ (*L,.(2;i, i')*L(a;3, t))H^ {n^{x2,t')-^\,,{xx,t')-^^Jx3,t)n^{xi, t)) ^^ 
and 

{B^{xi,X2,t')Bl{^X3,X4,,t)) u„ ~ 

U{x^,X2)U{x3, x^) U^l^^{xut')n^{x2,t')n^{x4, i)*L,.(a;3, t)) h^. (*R,a(a^i> ^')*L(a;3, t)) h^ 



+ (*[,.(a;i, i')*L,a(a^3, 0)ffL (^R.^Ca-'i, i')^R(a;2, i')"R(a;4, i)*L(a^3, 0)ffH 



(B14) 



C];{t',t) 

and 



(rcnL(xi, t')*R.a(a;2, t')*R,.(2^3, t)n^{xA, t))H^ 



Such that 



{h^{xx,t')'^^,^{x2,t')'^1 {X3,t)n^{xi, t))H^ 



L/+' +M< 



To evaluate the remaining three-particle correlations [C^(i',i)] 
it will be useful to introduce the following contour-order 
quantities 



- (*L a(^3, t)h^{xA, t)n^{xi,t')-^i^,^{x2,t'))H^ , 

(B18) 



and 



(B15) [Cf{t\t)\ 



nR(x-i, t')*R,^(2;2, t')*R.<T(a;3, t)h-p^{xi, t))H^ , 

(B19) 



(B16) [C^{t\t)]^^ 



- (*R,^(a;3, t)n^{xi, t)nR(xi, t')*R,<T(2;2, i'))if:,- 

(B20) 

The three-particle correlation functions (B15) and (B16) 
can be expressed as, neglecting two-particle correlations 



\^^') with three equal times and the three-particle vertex, 
I 



c^(i', t) « Y^^Tc^U, (.Ti, t')*L..' (X4, i))(rc*L,.' {xut')^^^A^2.t')^iAx3,t)ml^, [X,, t)} 

a' 

~ E(^C*L,<.(a^l, i')*[..(2^3, 0)(rc*[..(2:i, i')*L,.(2^2, i')*[,.' (2^4, O^L,.' (2:4, t)) 
a' 

+ E^^C*w(a^l,iO*L'(^4,t))(rcVl'L,(.Ti,t')*L..(^2,t')*L(a:3,i)*L..'(^4,0) 
a' 

- E^^C*L,.(^2, i')*U(^4, t)) (rc*[_,,(a;i, t')*L,.' (Xi, *')*[.. (a;3, t)*L,.(2;4, *)) 

cr' 

+ (rc'J'L,.(-^2, i')*[..(2;3, i))(rc^L(2;i, ^')^L(a^4, 0)- 



(B21) 



A similar expression exists for C^{t',t), Eq. (B16). 

Again assuming we are in the point contact or localized tunneling regime, where the wave functions of one side, say 
L, are spatially localized about the tunneling point, a;, the dominate contribution from the spatial integrals of (B2) 
can be approximated by the value at x. This is equivalent to neglecting the momentum dependance of the matrix 
elements of the interaction t/k,k' = (^L,k|t^|V'R,k}, which for a good metal is approximately true; especially for the 



10 



small energy region about the Fermi energy that one is typically interested in for tunneling experiments. Also for 
a contour-ordered quantity such as X{t,t') = Y{t,t')Z{t,t'), the greater and lesser functions are simply giving by 

[X(t,t')]^ = [Y{t,t')]^[Z(t,t')]'^. Along with defining 



/VCJ]^ ,CT"2 ,(73 ,£74 \ 



doj 



and the two-particle density of states pj/^(a;, w) = — ifmG]j^^(x,a;), where 
Eq. (B2) can finally be evaluated and is giving by 



[x,uj)e' 



-iu>(t-t') 



, (B22) 



(B23) 



eU^Yl J '^^duj' p^^^{x,UJ + eV)p^^^,{x,Lj'){[x^,^^^^^^,{x,uj' ~ uj ~ eV) + Xw',w.w.w'i^.^' - ^)]M^')[^ ~ M 
- [Xw,w',a',wix,^- w + eV) + x^-, -, -(x, w - io')]np{^ + eV) [l - np{^')] } 

X ^np{uj + eV)np{io') [l - np(a; + cj' + eV")] - [l - np{uj + eV)] [l - np{Lo')]np{uj + uj' + eV)\ 

-eirU^Y^ / dujduj' p^-^{x,uj + eV)p^,^{x,uj')p^^,^j^{x,uj+uj') 

X |np(a;-|-ey)np(w')[l -np(a; + a;')] - [l - np{Lo + eV)] [l - np(a;')]np(a; + u;')|, (B24) 



r 



which is Eq. (31). 
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